Coulomb interactions between disordered charge distributions 
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We review some of the recent results in the context of the Coulomb interaction between dielectric 
surfaces which are randomly charged but remain net-neutral on the average. Such surfaces are found 
to interact in vacuum with attractive long-range normal forces due to the presence of the quenched 
charge disorder which can compete with the standard Casimir-Van der Waals forces. They can also 
interact with random lateral forces and torques, which exhibit zero mean values but finite and long- 
range root-mean-square values. These effects can play an important role in Casimir experiments 
as well as in the interaction between solid surfaces as well as biomolecular surfaces that are often 
covered by disordered charge distributions. 



I. INTRODUCTION 



In most studies of electrostatic interactions between charged bodies a number of simplifying assumptions are made. 
Beyond the purely geometric simplifications, the charge distribution is often taken to be uniform. This assumption is 
clearly always an idealization as charge distributions in many systems will be inherently complex and/or disordered 
[IH5]. Examples of charge disorder are common in colloidal and soft matter systems [SHE]; specific examples include 
surfactant coated surfaces [3110] and random polyelectrolytes and polyampholytes [IT]. Metallic and dielectric surfaces 
with local dielectric constant variation can also exhibit charge disorder as local variations of the crystallographic axes 
of an exposed surface lead to a random surface potential, the patch effect [2HS1 . Finally the chemical preparation 
of samples is never perfect and charged impurities abound. The presence of charge disorder, even if the system is 
overall net neutral, can be shown to have strong effects on the interactions between bodies. Notably, charge disorder 
can lead to interactions which can mask the Casimir effect and may play an important role in the Casimir effect 
experiments, possibly making their interpretation rather delicate PHI]. 

The disorder we will consider in this chapter is defined as quenched as it is fixed once and for good in the preparation 
of the system, Fig. [TJ as in the case of charged impurities which are frozen in the boundaries of the materials and 
cannot move or react to electric fields acting upon them. It is for this reason that quenched disorder is sometimes 
referred to as frozen. In small systems, where the objects are not held fixed and can move and/or rotate with respect 
to one another, the system will tend to lower its electrostatic energy by aligning and positioning its components 
appropriately. This effect is believed to be a key component in the so-called lock and key mechanism, that plays a 
fundamental role in the biological recognition mechanisms [Ml [14] . 

In a typical Casimir force experiment where one has a two plate or a sphere-plate configuration, the surfaces in 
question are held fixed (not free to rotate or laterally to displace with respect to each other) and the two quenched 
charge distributions will be completely uncorrelated. If the sphere in a sphere-plate geometry is held close to one 
part of a large plate, the charge disorder will lead to random normal and lateral forces on the sphere, as well as 
random torques. If we carry out the same experiment between two different apposed parts of the sphere and the same 
surface, the fact that the sphere will feel a different charge distribution will lead to different normal and lateral forces 
and torques, in much the same way as if we had changed the sphere and plates by new ones produced by the same 
productions process, i.e., giving the same statistical disorder. Carrying out a sequence of such experiments will lead 
to an ensemble of measured forces, which can then be averaged to obtain the average force. In most cases we will 
see that the average normal force is non-zero. However, if the charge disorder in the plate is invariant by translation 
in space, that is to say it looks statistically the same everywhere across the plate, the average translational forces 
and torques will be zero. However for the translational forces and torques there will be a non-zero variance and they 
will fluctuate. The amplitude of these sample-to- sample fluctuations will give us an additional information about the 
nature of the charge disorder and may be useful in unravelling the various components of the force measured in typical 
Casimir or other force detection setups [2H1]- 
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FIG. 1: (Left) A schematic quenched surface charge distribution. (Right) Two parallel dielectric slabs Si and 52 which are 
semi-infinite in the z directions and infinite in the r = (x, y) plane are placed at a separation distance I. The disorder charge 
is taken to be distributed over the whole surface in slab S2 but is restricted to the finite region r £ A in slab Si. 



II. NORMAL ELECTROSTATIC FORCES BETWEEN CHARGE DISORDERED SLABS 

Consider a system of two parallel semi- infinite dielectrics Si and S2, with local dielectric constants e(x). We take 
e(x) = £1 in fiii, e(x) = £ 2 in S 2 and e(x) = e m in the intervening medium. We denote by I the separation between 
the two plates and by p(x) the quenched charge distribution for a given configuration of the dielectric bodies. The 
electrostatic energy is given by 

E(l) = \JJ dyrfxp(x)G(x,y;/)p(y) (1) 

where G(x,y; I) is the Green's function obeying 

e o V-[e(x)VG(x,y;0] = -«(x-y) ) (2) 

where £0 is the permittivity of vacuum. The Green's function depends explicitly on I as the overall spatial dielectric 
function depends on I. The average electrostatic energy of this configuration is thus given by 

(E(l)) = i J dxdyG(x, y; I) (p(x)p(y)) , (3) 

where (• • • ) indicates the disorder average over the random charge distributions. The total charge distribution is given 
by p(x) = pi(x) + P2(x), where pi(x) and p 2 (x) are the charge distributions in S± and S 2 . In some cases, notably 
for the estimation of lateral forces and torques, it is useful to assign one of the bodies, say Si, to have a charge 
distribution which is restricted to a subregion of Si which we will denote by A, see Fig. [T] 

Consider the case where charge distributions on different bodies are uncorrelated and where the average charge on 
each body is zero, i.e. (pi(x)) = (p2(x')) = and (pi(x)p2(x')) = for all x € Si and x' g S 2 . Then 

(E(l)) = l -J dxdyG(x,y;0[(pi(x)p 1 (y)) + (p 2 (x)p 2 (y))]. (4) 

We note that the average interaction has just two self interaction terms as there is no interaction on average between 
charges on different slabs, i.e., the interaction energy between a charge in Si with a randomly chosen charge in S 2 is 
on average zero as the charge in S 2 has an equal probability to have either the same sign or the opposite one as that 
in Si . A non-zero interaction energy is only possible between a charge in a slab and its image charge in the opposing 
slab, the sign of this image charge depending on the charge in question and the dielectric constants of the system. 

We notice, that there is an interaction between Si and S 2 , even if there is no charge in one of the slabs, stemming 
from an average interaction energy of a randomly charged slab (e.g., Si) with a charge- free slab (e.g., S 2 ). The origin 
of a net interaction between overall neutral surfaces is thus easy to understand, but still rather subtle and surprising. 

We now consider a simple model of the charge disorder, that can originate in either the surface charge or bulk 
charge. On the surface Si the surface charge distribution is 

Pi s (x)=^^(r-r i )^ + -) ) (5) 



i=l 



where r is the two dimensional in plane coordinate. If the area of the slab Si covered with charge disorder is A and 
the charges have the values qi = ±qi s with equal probability and N± s = n\ s A where n± s is the surface concentration 
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of charges, we easily find that (p s (r)) = and 

( Pls (x) Pls (x')) = g ls 5(z + l -)5(z' + l -)5(r r'), (6) 

where g\ s — nx s qf 8 . Typical values of q\ s are given by the electron charge e and the bulk impurity charge densities 
are between 10 10 to 10 15 e/cm 3 [TSJ[T5]. Using the size of a molecular layer we can thus estimate the typical surface 
charge density of the surface, generated by cutting the bulk. Experimentally the heterogeneous structure of the charge 
disorder on dielectric surfaces can be measured using Kelvin force microscopy pp. 
Away from the surface, we assume that bulk charge distribution has the form 

N u 

Pift( x ) = ^2<H$(x-Xi), (7) 

i=l 

where the charges qi take the value iqib with equal probability and have concentration m&. The bulk correlation 
function is then 

(p lb (x) Plb (x')) - gi b d(z - z')8{v - r'), (8) 

where gu, = nu,q\ b . The disorder assumed here has a short range correlation, an assumption easily modified to more 
general forms |17H19j . To give explicit formulas we will continue with short range correlation functions and assume 
that the (infinite) surface areas S 2 = S\ = S (by taking A to infinity in the case of slab Si) are completely covered 
by charge disorder. From Eq. Q we find that 

W)) = f [ffi.G(0, + g 2s G(0, 



S f — S r°° 

- dz g lb G(0,z,z;l) + - dz g 2b G(0, z, z;l) (9) 

1 J -00 1 Jl/2 

The Fourier transform in the in-plane coordinates is 

G(k, z, z'; l) = J dv G(r, z, z'; I) exp(-ik • r) (10) 

and we can use equation ^ to show that 

eo-^e(z) ^G(k, z, z'; I) - £ e(z)fc 2 G(k, z, z'; I) = -S(z - z') (11) 
dz dz 

where k = |k|. The above equation may be solved and we find that, for instance, when z < —1/2 we have 

= _J_ / [A^xp^O-A.exp^fcQlexp^A (12) 
v ' ' 2 £l e k\ 1- A 1 A 2 exp(-2/fcO J' y J 

where A a (for a = 1 or 2) is the dielectric jump parameter 

A a = ^-^. (13) 

With this result we can evaluate the terms appearing in equation Q using the inverse two dimensional Fourier 
transform 

G(0, z, z; = J dk G(k, z, z; 0- (14) 

For short range disorder, all the integrals can be evaluated analytically. It turns out that some integrals diverge due 
to self energies that are independent of the separation of the plates and therefore do not contribute to the force. 
The total average force can be decomposed into the contribution coming from each of the charge distributions in the 
system. For instance the average normal force (f^) due to the random surface charge on slab 1 is given by, 
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and the contribution of the force from the surface charge on slab 2 is obtained by switching the index 1 for 2 an vice 



a versa in Eq. ( 15 ) 



We first note that the interaction is an unscreened l/l 2 power law. Furthermore we see that the sign of the 
interaction is rather nontrivial and depends on the signs of the dielectric jump parameters A a [T71 HH1 [20] . We can 
make the following observations for, e.g., the force due to the slab 1, /f" . (i) for Ai > and A2 > 0, we find that 
< (attraction), (ii) for Ai > and A 2 < 0, (f^) > (repulsion) (iii) for A x < and A 2 > 0, < 

(attraction) and (iv) for Ai < and A 2 < 0, (f^) > (repulsion). So we see that the surface charge on slab 1 
feels an attractive force toward slab 2 if A 2 > 0, i.e., if the material composing slab 2 is more polarizable than the 
intervening medium. In the case of A 2 > but Ai < 0, the force due to the charge on slab 1 is attractive, but due 
to the charge on slab 2 is repulsive! It is interesting to note that the sign of the thermal van der Waals interaction 
between two dielectric slabs depends only on the product AiA 2 as we shall see below. 

The normal force due to the bulk charge disorder in slab 1 can be shown [T7] to be given by 



(/<?>> 5l6 A 2 (l-A?) 



167T£i£ Z 1 - AiA 2 



(16) 



The result for the force due to the charge distribution in the slab S2 is obtained on switching indices. The sign of 
the force depends purely on that of A 2 (Ai) andbehaves in the same way as the surface charge. However, we see 
that the force generated by bulk charge disorder is longer ranged than due to the surface charges, decaying as l/l. 
In the above computation it is important to bear in mind that slabs will have a finite thickness-changing the inter 
slab separation should not introduce/remove quenched charge from the system, which will give an effective bulk term 
in the energy. A way of avoiding this, without complicating the dielectric problem, is to assume that the quenched 
charge disorder has a support which is in the interval \—{L + 0/2, — J/2]> m slab 1 for instance, while keeping L fixed 
upon changing I. 

The components of the force due to both bulk and surface charge disorder in the same slab have the same sign, 
this means that the electrostatic component of the normal force is monotonic, having the same sign for all values of 
I. If one takes into account thermal Casimir forces which are also monotonically attractive or repulsive, it is possible 
to find systems where the overall force becomes zero, corresponding to a stable or unstable equilibrium jTHJ [TS] . The 
normal thermal Casimir or the zero Matsubara frequency van der Waals (vdW) force interaction behaves as 



(17) 



/vdw _ fc B TLi 3 (A 1 A 2 



S 8nl 3 

where Lis(z) is the trilogarithm function. If the prefactor is chosen to be negative (by taking AiA 2 < 0) the repulsive 
vdW interaction can stabilize the interaction, preventing collapse at short distances. It is theoretically possible to 
have an attractive charge disorder interaction, behaving as l/l or l/l 2 at large separations, meaning that the overall 
interaction potential will possess a stable minimum [T51 Qj|] . 

The above results can be easily generalized to the case where the charge disorder possesses a finite correlation length 
in the plane of the slabs [TTJ [TH]. However the case of a non-zero correlation normal to the slabs, in the z direction, 
is more difficult to deal with. 

Finally we should note that the physical usefulness of an average value depends on the scale of fluctuations of the 
force with respect to the average force. For simplicity we will consider the variance of the normal force with just 
surface charge disorder. After a rather long computation for surface charge disorder with zero-correlation length, we 
find the variance of the total interaction force induced by surface charge disorder to be given by [20] 



<(/i R) ) )c = 7-272 /2 (alDu + gl s D22 + 2g ls9 2 S D 21 ) (18) 



where 



A,- A; A 



£ 4 



« 1 3 -I-T.V2 



3(e m +£i) 2 (e m + £j) 2 



2A : A 
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(1-AiAa) 2 v 



(19) 



where we have define otij = (1 + <%)(— 1)* +J for i,j = 1,2 (note that no summation over the indices i,j is assumed in 



Eq. (19)). The normal force fluctuations due to surface charge disorder then scale as 



a/w „ il. (20) 
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The ratio of the force fluctuations to the average force R = Sfl , / (/i ) thus scales as R ~ l/y/S and are small as the 
surface area is increased (as expected from self averaging in the thermodynamics limit), however they increase as the 
inter-slab separation I is increased. In most Casimir type experimental setups the linear dimensions of the interacting 
objects are typically much greater than their separation, keeping the force fluctuations small [3J. 



III. LATERAL ELECTROSTATIC FORCES BETWEEN CHARGE DISORDERED SLABS 

Suppose now that the region containing the charge disorder in the dielectric slab 1 is laterally smaller than the slab 
2 whilst remaining large of area A, see Fig. [I] We furthermore assume that charge disorder is distributed only on the 
bounding surfaces of the two slabs. If we displace slab 1 laterally by a vector a within the plane of the slabs, this will 
induced a change in the charge distribution 5p(x) and give a change in energy 

6E = J dxdy6p(x)G(x,y;l)p(y). (21) 

Note that in this set up we assume that there is no change in the dielectric function and so the only displaced charges 
are the quenched ones in slab 1. If we move slab 1 by an infinitesimal amount a it is clear that the only non-zero 
contribution to Sp is 8p\. The charge density associated with the slab 1 can be assumed to consists of point charges 
qi at the points Xj = (r^, Zi) within the slab 1 (or on its surface in the case or surface charges). If the change in the 
system is to displace slab 1 by an infinitesimal vector a, the new charge distribution associated with slab 1 corresponds 
to all the charges at being displaced to (r^ + a, zi) when the whole body is displaced by a. This then gives 

<5p(x) = -a-V r pi(x). (22) 

The change in energy on displacing body i through a is given only by the interaction of the charges and image charges 
in slab 1 with those in 2, and therefore defines a lateral electrostatic force f that allows us to write the average of this 
interaction energy as 

(<5£(a)) = -(f«.a) = -a. j dxdy (V Pl (x)G(x, y; l)p 2 (y)), (23) 

The average lateral force is thus zero as we assume that charge distributions on different bodies are uncorrelated. 

Nevertheless the variance of the lateral force will be non-zero and can be computed using the two point correlation 
function of the charge disorder. For brevity we state only results for surface charge disorder with spatially uncorrelated 
charge disorder [20,. We find that the lateral force components in the in-plane directions i, j = 1,2 have the correlation 
function 

/f(0 f (0\_ ASjjg^sel, , AA^ (OA\ 

ifi f * } ~ -4«e$P(e m + + ea )=»A 1 A 2 H ' " A ^ (24) 

showing that the typical lateral force fluctuations scale as [30] 

f/ W ~ ^, (25) 

the same scaling as for the fluctuations of the normal force. We note that in order for the variance to be non-zero we 
do not require any dielectric discontinuities and find that for s\ — £2 — £m 

I M) Al\ _ M ij9u92s , 0( ,, 



IV. ELECTROSTATIC TORQUES BETWEEN CHARGE DISORDERED SLABS 

On rotating the slab 1 by an angle around its symmetry axis [21] . that is to say in the direction perpendicular 
to the normal between the bounding surfaces of the two dielectric media, we find (assuming that charge disorder is 
distributed only on the bounding surfaces) that the new charge distribution of slab I is given by 



p'i(x)= ^2 q n S{r - R e r n )S(z ~ z n ), 

neS! 



(27) 
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where Rg is the two-dimensional rotation matrix. For an infinitesimal rotation angle 69, one has Rse = 1 — i(S9) 02, 
where a% = ( ' ] is the Pauli matrix. The only non-zero contribution to Sp in this case (assuming the summation 



x 1 , 

over the in-plane Cartesian components i,j — 1, 2) is then 



5p(x) = t(66) («7 2 ) u -r,— Pl (x). (28) 



so that the change of the interaction energy, Eq. (21 ), on rotating the surface Si, is due to the interaction of charges 



and image charges in Si with those in S2. We may thus write 

SE = t(56)(a 2 )ij I dr'drdzdz'lr'j—p^r'sz')] 

xG(v-v',z,z'-l)p 2 {v,z), (29) 

where r' and r are again the two-dimensional coordinates in the planes of Si and S2 respectively and z' and z' are 
the coordinates normal to the planes. We again assume that the charge disorder is restricted to a subregion of area 
A of Si, see Fig. [I] In an experimental set up, it would be more likely that Si is finite and the charge distribution is 
over all Si, however, in the approach followed here we avoid dielectric edge effects and the Green's functions involved 
then retain their planar symmetry. 

Integration over the coordinate r' is over the finite area A, while that over r is unrestricted. The torque r acting on 
the surface Si is then given by SE = —{86) r. As the charge distribution on the surfaces Si and Si are uncorrelated 
we find that (SE) — — (S6)(t) — 0. Thus the mean torque is zero. 

From here one can derive |21j an intuitively clear physical relation between the lateral force fluctuations and the 
torque fluctuations of the form 



r 



2 > = j(ti l) fi l) )Iii, (30) 



where the moment of inertia tensor is defined as 



Iij = I dr (Sijr 2 - nrj). 



A 



The torque fluctuations are thus connected with the lateral force fluctuations through a geometric factor encoded by 
moment of inertia tensor of the region A; thus the torque fluctuations are not the same for a square or a disc of the 
same surface area. This result is completely general and valid for the assumed plane-parallel arrangement of the two 
disorder-carrying dielectric surfaces. 



As the lateral force fluctuations have the scaling y A/l, (25 1 , we find that 



(if /f ) ~ ~. (31) 



It then follows that typical torque fluctuations scale as 

8t~j, (32) 

being extensive in terms of the area A. This is as one might expect because torque is determined by the geometry 
of the area A even in the limit of large area (which is not the case for the random lateral force derived above) . The 
geometry-dependence of the torque fluctuations and the scaling with area A are obtained simply from the moment of 
inertia tensor. 

Again we note that in order for the torque variance to be non-zero we do not require any dielectric discontinuities 
and find that for ei = £2 = £m the torque fluctuations are given by 

/ 2\ _ 9ls92sA 2 

{T } ~ I28^ 2 e%e%p ' ^ 
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V. LESSONS 



Charge disorder has important repercussions for electrostatic interactions between net-neutral dielectric slabs bear- 
ing random charges on their bounding surfaces and/or in the bulk. Quenched disorder leads to an additive contribution 
to the net (normal) interaction force that scales as l/l (or I /I 2 ) for bulk (or surface) charge disorder, may be attractive 
or repulsive and depends on the dielectric contrast of the materials. 

Because of the nature of electrostatic interactions between disordered media, there is a sample-to-sample variance 
in the lateral as well as normal forces and it can be substantial. In the ideal limit where the probe area is very large, 
the average force is much larger than its fluctuations. However, in some experimental setups the probe size may 
be quite small and the sample-to-sample force fluctuations could become important. In addition, force fluctuations 
become important at large separations where the normal force is weak. In the case of lateral force variance, since 
the average is zero, fluctuations are the only thing observable. Fluctuations in the normal and lateral directions are 
always comparable; for the special case of a uniform dielectric constant the variance of the force fluctuations in the 
normal direction is exactly twice the magnitude of the one in the lateral direction |20) . 

The disorder-induced torque fluctuations scale differently with the surface area of the interfaces bearing charge 
disorder, with the typical torque fluctuation being extensive in the surface area A. This opens up a feasible way to 
measure charge disorder-induced interactions between randomly charged media, in a way which is independent of 
normal force measurements and with a higher signal-to-noise ratio than lateral force measurements. 
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